Tensor effective interaction in self-consistent Random Phase Approximation 

calculations 



M. Anguiano\ G. Co' 2'^, V. De Donno^'^ and A. M. Lallena^ 

1 ) Departamento de Fisica Atomica, Molecular y Nuclear, 
Universidad de Granada, E-18071 Granada, SPAIN 
2) Dipartimento di Fisica, Universitd del Sa lento. Via Arnesano, 1-73100 Lecce, ITALY 
3) INFN, Sezione di Lecce, Via Arnesano, 1-73100 Lecce, ITALY 
O ! (Dated: January 25, 2013) 



> 



We present a study of the effects of the tensor-isospin term of the effective interaction in Hartree- 
Fock and Random Phase Approximation calculations. We used finite-range forces of Gogny type, 
and we added to them a tensor-isospin term which behaves, at large internucleonic distances, as the 
analogous term of the microscopic interactions. The strength of this tensor force has been chosen 
to reproduce the experimental energy of the lowest 0~ excited state in ^®0, which shows large 
\l ' sensitivity to this term of the interaction. With these finite-range interactions, we have studied the 

effects of the tensor-isospin force in ground and excited states of carbon, oxygen, calcium, nickel, 
zirconium, tin and lead isotopes. Our results show that the tensor force affects mainly the nucleon 
single particle energies. However, we found some interesting cases where also bulk nuclear properties 
are sensitive to the tensor interaction. 

Qji PACS numbers: 

I. INTRODUCTION 

K*" ■ At the beginning of the 40's of the past century, the existence of the electric quadrupole moment of the deuteron 
I [II) [2| was explained by Rarita and Schwinger by introducing a static tensor term in the nuclcon-nuclcon (N-N) 
interaction Since then, tensor terms are unavoidable ingredients of the microscopic N-N interactions, i.e. those 

interactions constructed to reproduce the properties of two-nucleon systems. 

Despite their relevance in microscopic interactions, the tensor terms are usually neglected when effective interactions 
and theories are used. In the effective theories some complicated many-body effects are treated, obviously in effective 
and average manner, by changing the values of the parameters of the interaction. Some specific observables of the 
nucleus are chosen to select these values. For example, in Hartree-Fock calculations these observables are usually 
the nuclear binding energies. The effective theory is expected to be able to describe other observables. If this fails, 
one searches for many-body effects which should be explicitly treated to improve the description of the data. In this 
manner, we link many-body effects to specific observables. In the case of our interest here, the tensor force, the point 
5^ is the identification of observables clearly depending on the presence of this term in the effective N-N interaction. 

In these last years, the interest on the tensor terms of the effective N-N interaction has increased because the 
inclusion of these terms improves the description of the single particle (s.p.) energies of some isotope or isotone chains 
when the Hartree-Fock theory is used 0-[l5| . 
We see some weak points in using the s.p. energies to define the strength of the tensor terms of the effective 
interactions. First, s.p. energies are extremely sensitive to the spin-orbit terms of the N-N interaction, and this 
obscures the possibility of a clear identification of the tensor effects (see for example the discussion in Ref. [l6|). 
Furthermore, observations are always done on global nuclear properties, therefore the identification of the measured 
quantities with s.p. properties of the nucleus is done by imposing to the observed quantity the physical interpretation 
given within a mean field description of the many-body system. The fact that experimental values of spectroscopic 
factors are usually rather different from the mean field expectations is a clear indication of the limits of this procedure. 

In this article we propose an alternative approach to select the strength of the effective tensor forces. We have 
looked at the excitation spectrum to find observables particularly sensitive to the tensor force. We have identified 
these observables with the energies of the 0~ charge conserving excitations. The large sensitivity of the 0~ excitation 
to the tensor parts of the nucleon- nucleon interaction was pointed out already in Ref. (l7j . By using a recursive 
self-consistent Hartree-Fock (HF) plus Random Phase Approximation (RPA) procedure, we have chosen the strength 
of the tensor term of the effective interaction to reproduce the experimental value of the 0~ in ^^O. With these new 
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interactions we have investigated the ground and excited state properties of various nuclei by doing HF plus RPA 
calculations. 

The structure of the tensor terms we have considered and the methodology used to select their strengths are 
presented in Sec. |TT1 We discuss the results obtained in the description of the ground states of various nuclei in Sec. 
Illli and in Sec. |TV]the results obtained for the excited states. We summarize the main points of our work and draw 
our conclusions in Sec. |Vl 



II. THE INTERACTION 



The most important tensor component of the microscopic N-N interaction is that related to the tensor-isospin 
channel [3, [l^ whose long range behaviour is dominated by the exchange of a single pion. Since the pion is the 
lightest meson, the range of the tensor-isospin term is longer than the ranges of the other terms of the N-N interaction. 
For this reason we have chosen to consider, in our effective interactions, only tensor-isospin terms with finite range. 
The use of finite range forces requires, in both HF and RPA calculations, the evaluation of direct and exchange 
interaction matrix elements. 

The tensor-isospin term of our effective interaction is based on the analogous term of the microscopic Argonne 
V18 interaction [1^. We have multiplied the radial part of this term by a function which simulates the effect of the 
short-range correlations [l^l . In our work, the radial part of the tensor-isospin term has the form 

ve{r) = ve.AYisir) [l - exp(-6r^)] , (1) 

where we have indicated with r the distance between the two interacting nucleons, with fg.Avis the radial function of 
the Argonne V18 tensor-isospin potential fiot . and with b a free parameter. The changes in the tensor-isospin term 
produced by choosing different values of b are shown in panel (a) of Fig. [TJ where we present the Fourier transformed 
function Ve{q), defined by the equation 

Vg (q) Si2 (q) = y" d V exp (z q • r) (r) ^12 (r) 

= -Att drr'^j2{qr)ve{r)Si2{r). (2) 



In the above equation r and q indicate the moduli of r and q, and we used the definition of the tensor operator 

5„M^3 '-'"-'-lp'--l -.(l)-.(2). (3) 

where cr are the usual Pauli spin matrices. 

The results presented in the panel (a) of Fig. [T] show that the correlation effects reduce the strength of the bare 
tensor force. The smaller is the value of b, the more extended in r space is this effect, therefore the function to be 
integrated in Eq. ([2]) becomes smaller. 

The first step of our study consisted in identifying an observable very sensitive to the tensor force with the aim 
of using it to determine the strength of this part of the effective interaction. We have conducted this study by 
investigating the excitation spectra of the ^^C, ""^^O, '^"Ca, ''®Ca, ®°Zr and ^°*Pb nuclei within the phenomenological 
RPA approach developed by the Jiilich group [2l|, [23. In this approach, based on the Landau-Migdal theory of finite 
Fermi systems [l^, the set of s.p. energies and states is obtained by using phenomenological mean field potentials 
that, in our case, have the shape of Woods-Saxon wells. The values of the parameters of the potential are chosen for 
each nucleus in order to reproduce at best the empirical values of the s.p. energies around the Fermi surface and the 
values of the charge root mean square radii. The explicit expression of the potential, and the values of the parameters, 
can be found in Refs. [l^, [13]. Following the philosophy of the Landau-Migdal approach, in the RPA calculations we 
substituted the Woods-Saxon s.p. energies with their experimental values, when they are available. 

All the RPA results presented in this article have been obtained by using a discrete s.p. basis. In the phenomenolog- 
ical calculations discussed in this section we have used a discrete s.p. basis obtained by diagonalizing the Woods-Saxon 
well in a harmonic oscillator basis. In analogy with the work of Ref. [13], we have used configuration spaces large 
enough that the inclusion of additional s.p. states does not modify the energies of the first excited states, below 20 
MeV in the lighter nuclei and 15 MeV in the heavier ones, within 0.1 MeV. 

In these phenomenological calculations we used as basic N-N effective interaction a density dependent Landau- 
Migdal force. The values of the force parameters, given in Refs. 0, HI], change for each nucleus, and are chosen to 
reproduce at best the energies of the collective low-lying 3~ states of each nuclei, and also the energies of the 12~ 
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Figure 1: Momentum dependent term of the tensor force, see Eq. ([2]), for various parameterizations. In the panel (a) we 
compare the tensor-isospin term of the bare Argonne V18 interaction (full line) with the interactions obtained from Eq. ^ 
by using different values of the parameter b. In panel (b) we compare the bare Argonne V18 term with the tensor term of 
the GT2 force of Ref. O (dotted line) and with that of the M3YP2 force of Ref. [38] (short-dashed line). The tensor 
forces we have constructed, named DIST and DIMT, are indicated by the dashed-dotted and dashed-doubly-dotted lines. The 
long-dashed line (CBF) shows the tensor term obtained by multiplying the bare Argonne V18 interaction with the scalar part 
of the correlation function obtained in microscopic Correlated Basis Function calculations [20l ]. 





excited states in ^"^Pb. We added to this basic interaction a tensor term of the form given by Eq. ([T]), and we studied 
the excitation spectra of the nuclei mentioned above. 

In agreement with previous calculations [22, 26] we observed that the influence of the tensor term on the natural 
parity states is negligible. We found few cases of interest in the spectrum of the unnatural parity excitations. The 
most interesting one was the excitation of the 0~ states, which show common characteristics in all the nuclei we 
considered. 

In Fig. [5] we have summarized the results obtained for the energies of the lowest 0~ excitations in the nuclei 
investigated. In this figure, the squares indicate the values of the 0~ excitation energies, obtained by adding the 
tensor term of Eq. ([Ij to the Landau-Migdal force, as a function of the parameter b. The horizontal dashed lines 
indicate the values obtained without tensor term, while the dotted lines show the values obtained by using the full 
tensor term of the Argonne VI 8 interaction. The full lines show the experimental values (27l. [28j. 

We observe that the effect of the tensor term is always attractive, i.e. all the energies obtained with the tensor 
terms are smaller than those obtained without it. The values of these energies decrease, monotonically and smoothly, 
from the dashed to the dotted lines when the value of b increases, as we have naively expected. Discussing the results 
in more detail, we have obtained variations of the energies from 4 up to 8 MeV, and, in relative variations, between 
30% and 100%. These are the largest effects of the tensor force on the excitation energies we have found in our 
investigation. In the ''^Ca case, the effect of the full tensor term is so strong that we obtained an imaginary solution 
for the RPA equations. 

The remarkable sensitivity of the 0~ excitation energies to the presence of the tensor force, and their smooth 
behaviour with respect to the changes of its strength, make these energy values particularly suitable to be chosen as 
experimental benchmarks to select the strength of the tensor terms in effective interactions. 

In this work we have performed HF plus RPA calculations. The HE equations were solved with the method used in 
Refs. [l^,!!^]- This method is based on the plane wave expansion technique developed by Guardiola and Ross [3ll.[3^. 
After the iterative process has reached the convergence, we solved the HE differential equations not only for the s.p. 
states below the Eermi level, the hole states, but also for those states above the Eermi level, the particle states. The 
numerical method automatically produces a set of discrete levels even when the s.p. energies are larger than zero, i.e. 
in the continuum region. We did not find general criteria for the stability of our results. This problem is not very 
important in the phenomenological RPA approach, since the effects related to the truncation of the s.p. configuration 



4 



20.0 
IS.O 
16.0 
14.0 
12.0 

12.0 
TO O 

3 

8.0 



' I 


9.0 








7.0 




" ■ ■ 


5.0 






J.O 







7.0 














5.0 










- '"'Zr ■ ■ ■ .- 






3.0 





;oo 




6.0 




s.o 


:" ■ -'"ca : 


4.0 




6.0 


" ■ ■ . 










2.0 




4.0 









0.0 0.4 as 1.2 1.6 0.0 0.4 OS 

b [fm-2] b [fir 



Figure 2: The squares show the values of the excitation energies of the lowest 0~ states for various nuclei as a function of 
the value of the parameter b of Eq. ((TJ ruling the strength of the tensor force. The full lines show the experimental values, 
the dashed lines the values obtained without tensor, and the dotted lines the values obtained with the full tensor term of the 
Argonne V18 interaction. 



space are taken into account by changing the parameters of the force. However, in self-consistent calculations this 
is a more serious problem, since the interaction parameters are chosen to reproduce, in HF calculations, the ground 
state properties of the various nuclei, and the force remains unchanged in RPA calculations. 

To keep under control these problems, we restricted our study to the low-lying excited states. For each excitation 
studied we used the same criterion considered in the phenomenological approach, i.e. we controlled that the energy 
eigenvalues of the first low- lying states did not vary by more than 0.1 MeV against the enlargement of the configuration 
space. In order to obtain this numerical accuracy, we had to use configuration spaces composed by a few thousands s.p. 
states, more than 2000 in ^°*Pb. The numerical stability of higher energy excitations, such as giant resonances, requires 
even larger configuration spaces. In this case we believe it is necessary to abandon the discrete RPA calculations and 
treat correctly the continuum, as it is done, for examplej_in Refs. 33, 34]. 

We have built two new forces by adding to the DIS [33] and DIM ^] parameterizations of the Gogny interaction 
[37} a tensor- isospin term similar to that given by Eq. ([T]). We label DIST and DIMT these new interactions. Since 
in Gogny-like forces the spin-orbit term is fixed to reproduce the experimental splitting of the s.p. energies of the 
lpi/2 E^nd lp3/2 neutron states in ^^O, we used this nucleus as reference. The new interactions have been fixed by 
using an iterative procedure. We started with a HF calculation without tensor force to produce a set of s.p. energies 
and wavefunctions to be used in the RPA calculations. Then, we made a RPA calculation with the tensor force and 
we fixed the value of the parameter b of Eq. ([1]) in order to reproduce the energy of the first 0~ excitation of ^^O at 
10.6 MeV. With this new interaction we recalculated the HF s.p. energies by changing the spin-orbit interaction to 
reproduce the splitting of the two p states quoted above. These HF and RPA calculations have been repeated until 
the convergence of the result was obtained. By using this procedure we have found for the parameter b the value 
0.6 fm~^, for the DIST force, and of 0.25 fm~^, for the DIMT one. Summarizing, we added a tensor term to the DIS 
and DIM Gogny-like interactions and we modified only the spin-orbit terms from 130 MeV, in the original DIS force, 
to 134 MeV, in the DIST interaction, and from 115 MeV, in the DIM force, to 122.5 MeV, in the DIMT interaction. 
No other values of the force parameters have been changed. 

The tensor terms of the DIST and DIMT interactions are indicated in panel (b) of Fig. [l]by the dashed-dotted 
and dashed-doubly-dotted lines respectively. In this figure they are compared with the tensor-isospin term of the 
microscopic Argonne V18 interaction (solid curve). By construction, the effective tensor terms are smaller than that 
of the bare N-N interaction. More interesting is the comparison with the long dashed line which has been produced 
by multiplying the bare interaction with the scalar part of the short-range correlation function obtained in Correlated 
Basis Function calculations [l^. The remarkable difference between this line and those of the DIST and DIMT forces 
indicates that our procedure includes in the effective tensor term not only the effects of the short-range correlations, 
but also some other many-body effects that the microscopic calculations consider explicitly. In the same panel we 
make a comparison with other two tensor terms of finite-range interactions used in the literature, the GT2 [lo| and 
the M3YP2 [s^ forces. The tensor term of the GT2 is constructed to have the same volume integral of the Argonne 
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Figure 3: Panel(a): relative percentage difTerences, Eq. between binding energies calculated with the DIST and DIS forces 
(open squares) and with the DIMT and DIM (solid circles) interactions. Panels (b) and (c): binding energies per nucleon 
calculated with the various interactions compared with the experimental values (solid triangles) taken from Ref. [13, |4l[. The 
experimental values for the ^*0, ''°Ca, **Ni, ^®Ni and ^""Sn nuclei are estimated [IJ]. The lines have been drawn to guide the 
eyes. 



V18 tensor force. The strength of this tensor force is much larger than those of the tensor forces we have built. On 
the opposite, we observe that the tensor term of the M3YP2 force is much smaller. 

Even though the DIST and DIMT forces reproduce the experimental value of the excitation energy of the 0~ state 
at 10.6 MeV in ^^O, they produce rather different RPA wave functions, as we have verified by calculating transition 
densities and inclusive neutrinos cross sections for this excited state. 



III. HARTREE-FOCK RESULTS 



We used the DIST and DIMT interactions, whose construction has been described in the previous section, to make 
spherical HF calculations for a set of nuclei in different regions of the isotope chart. We have chosen nuclei where the 
s.p. states below the Fermi level are fully occupied to avoid deformation problems and to minimize pairing effects. We 
have verified these features by controlling the results given by the deformed Hartree-Fock-Bogolioubov calculations of 
Ref. [39]. 

In Fig. [3] we compare the binding energies obtained for various nuclei we have considered, with the experimental 
values taken from Refs. [l^, l4ll. The experimental energies of the ^^O, '^^Ni, ^''Ca, ''®Ni and ^"''Sn nuclei have not 
been measured but estimated |4ll |. 

The methodology used in our study is already evident. We have investigated the effects of the tensor force by 
comparing the results obtained by using forces with and without tensor term. To emphasize the effects of the tensor 
force, we show in the panel (a) of Fig. [3] the quantity 

Ag^lOO ^°^"^~^°^" , (4) 

which is the relative percentage differences between the binding energies calculated by using interactions with and 
without tensor term, SdiqT and Sdiq respectively (a = S,M). In this panel the open squares indicate the results 
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Figure 4: Relative percentage difference, Eq. ([5}, between root mean square radii of neutrons (upper panel) and protons (lower 
panel) distributions calculated with and without tensor for DIST and DIS (open squares) and DIMT and DIM (solid circles) 
interactions. 



obtained with the DIST and DIS interactions and the solid circles those obtained with the DIMT and DIM ones. 
The lines have been drawn to guide the eyes. 

If we exclude the anomalous values obtained for the ^''O nucleus, we observe that all the other results lie in the 
range AE « ±2. In general, the inclusion of the tensor term in the DIM force produces more binding, while it has 
opposite effect in the DIS case. These results show that effect of the tensor on the binding energy of the nuclei we 
have investigated is rather small, confirming the results of Ref. [s^. As a consequence, and as we can see in panels 
(b) and (c), the inclusion of the tensor term does not modify the agreement with experimental data in a significant 
manner. 

We have investigated the effect of the tensor force on the proton and neutron density distributions obtained from 
our HF calculations, and, also in this case, we found small effects. We summarize the results of this study in Fig. 
3] where we have shown the relative percentage difference between root mean square (rms) radii obtained by using 
interactions with and without tensor force, 

Ar = 100 ^^'''^"^ "^ ' ^ ^'"^^ . (5) 

We show the results for the neutron radii in the upper panel of the figure, and those for the proton radii in the 
lower panel. We have indicated with the open squares the results obtained with the DIS forces and with the solid 
circles those obtained with the DIM ones. 

The effects of the tensor force are rather small, of the order of few parts on a thousand. We observe again the 
different sign between DIS and DIM results. The larger binding produced by the tensor in the DIM case generates 
more compact nuclei, i.e. with smaller rms radii. In the case of the DIS forces the effect is just the opposite. This 
trend is present in both neutron and proton cases. 

The results we have presented so far indicate that the bulk properties of the nuclear ground states are not greatly 
affected by the presence of the tensor force. The situation changes when the s.p. energies are considered. In the 
remaining part of this section we shall discuss results concerning s.p. properties. Henceforth we shall distinguish the 
proton and neutron s.p. levels by using the tt and v labels, respectively. 

A first quantity we have studied is the difference between the s.p. energies of spin-orbit partner levels 

S = £(-1/2 - £(+1/2 ■ (6) 

In particular, we have studied the difference between the values of s obtained by using forces with and without tensor 
terms 

As = sdiqT - SDiQ ■ (7) 

In Figs. [5] and [S] we show the values of As calculated for the Ip, Id and 1/ proton and neutron levels, respectively, 
for all the nuclei considered. In these figures the open squares indicate the results obtained with the DIS interactions. 
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Figure 5: Differences between the s.p. energy difTerences of spin-orbit partners levels, Eq. ((7]), calculated with and without 
tensor forces, for the Ip, panel (a). Id, panel (b), and If, panel (c), proton states. The results for DIS and DIST interactions 
are shown by open squares and solid circles, respectively. The arrows indicate those isotopes where the effect of the tensor is 
expected to be zero. 



and the solid circles those obtained with the DIM ones. The arrows indicate those nuclei where all the spin-orbit 
partner levels, for both protons and neutrons, are fully occupied. 

Let us consider first Fig. [Sj where, for each isotope chain we show the evolution of As values with the increasing 
number of neutrons. We first observe that, for each s.p. level investigated, the DIS and DIM results have identical 
behaviour. Minima and maxima are in the same position for both type of calculations. A second observation is that, 
in general, As is negative. This means that the tensor force reduces the energy difference between spin-orbit partner 
levels. A third observation is that the effects of the tensor force are minimal, almost zero, for those nuclei indicated 
with the arrows. 

The first observation indicates that the effects we have pointed out are strictly related to peculiarities of the tensor 
force and of the nucleus investigated. The small difference between the DIS and DIM results reflects the difference 
between the tensor forces in DIST and DIMT, as we have shown in Fig. [TJ 

The second and third observations are well understood within the scheme proposed by Otsuka et al. [l^, 112] ■ 
The effect of the tensor interaction between a proton and a neutron occupied s.p. levels is attractive if one of the 
levels has an angular momentum j> =1 + 1/2 and the other one = I — 1/2. If the angular momenta of the two 
s.p. levels are of the same type, i.e. both j> or both j<, the effect of the tensor force has opposite sign. The effect 
of an occupied neutron level with angular momentum j> is to increase the energies of the s.p. proton levels with j>, 
and to lower those of the levels with j< . As a consequence, the splitting between the energies of the proton spin-orbit 
partner levels is reduced. This effect is reversed when the occupied neutron level has j<. If both j> and j< neutron 
levels are occupied, the two effects cancel with each other. 

The results presented in Fig. [5] are well explained within this picture. In the nuclei marked with the arrows, all the 
j> and j< neutron levels are occupied. In these nuclei we do not expect any tensor effect. In reality the values of As 
are not exactly zero even in these cases because we have changed the strengths of the spin-orbit interactions in the 
forces with tensor terms (see Sec. |ll]). In all the other nuclei there is, at least, one occupied neutron level with j>, 
whose spin-orbit partner level is empty. The effect we have discussed above predicts negative values of As, as those 
shown in the figure. 

In Fig. [5] we have shown the values of As for the neutron Ip, Id and 1/ levels. Since in the oxygen and calcium 
isotopes all the proton s.p. partner levels below the Fermi surface are completely occupied, the effects we observe in 
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Figure 6: The same as Fig. [5] for neutron states. 



these isotopes are due to the tensor interaction acting within neutron s.p. levels only. This is exactly the effect we 
have discussed above, but acting between s.p. states with the same isospin. This effect is weaker than that between 
states with opposite isospin, as the results for the ^^O, ^"^O, ''^Ca and ^■^Ca nuclei indicate. In these nuclei one of 
the neutron spin-orbit partner levels is unoccupied and the tensor effect is expected to be present. Actually, there is 
also an effect produced by the different values of the spin-orbit forces, but we have verified that, for the cases under 
consideration, this effect is negligible. 

The analysis of the As results for the Ni, Zr, Sn and Pb isotopes is more complicated, since, in these cases, in both 
the proton and neutron sectors, there are spin-orbit partner s.p. levels not fully occupied. For this reason, we must 
consider the interaction with both proton and neutron levels. The similarity between the behaviour of the results 
obtained with the different interactions indicate that we are observing an effect due to the peculiarities of the tensor 
force. 

Another quantity of interest related to the s.p. energies is the gap, g, between the energies of the levels just above 
and just below the Fermi surface. In Fig. [7] we have shown the proton (lower panels) and neutron (upper panels) 
energy gaps calculated for the nuclei we have investigated with the four forces we are using. Our results are compared 
with the experimental values (solid triangles) extracted from the binding energies of nuclei with atomic numbers 
differing by one unit. 

The first remark related to the results shown in Fig. [7]is that, in general, the effects of the tensor force are relatively 
small and they are similar for both type of interactions. They are negligible for neutrons, while some noticeable effects 
are present in the case of protons. The neutron results, obtained with and without tensor, follow reasonably well the 
behaviour of the experimental energy gaps. 

The effects of the tensor force are better emphasized in Fig. |S] where we show the difference 
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between the gap values obtained by using interactions with and without the tensor force. Also in this case we have 
indicated with the arrows the nuclei where all the spin-orbit partner levels are fully occupied, for both protons and 
neutrons. The results obtained with the DIS interactions are shown by the open squares and those obtained with 
the DIM interactions by the solid circles. It is interesting to observe that also in this case the behaviour of the DIS 
and DIM results is similar, though the effects produced by the DIST force are slightly larger than those obtained for 
DIMT. 

Also the results in Fig. [5] can be well explained within the scheme proposed by Otsuka. In the spin unsaturated 
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Figure 7: Neutron, panels (a) and (b), and proton, panels (c) and (d), energy gaps, in MeV, for the various nuclei and 
interactions we have investigated. In the left panels we show the results obtained by using the DIS interactions and in the right 
panels those obtained with the DIM forces. The experimental values, solid triangles, have been extracted from the binding 
energies of the neighbouring nuclei [i^, |4l[ . 



oxygen isotopes, the unpaired neutron levels are always of j> type. The tensor force lowers the energy of the {lpi/2)-n 
occupied proton level (j< type) and increases that of the (lc?5/2)ir empty level (j> type). For this reason, the results 
for the oxygen isotopes, in the proton case, have positive values. An analogous effect is present also for the "^^Ca and 
^^Ca isotopes. In this case, the states to be considered are the holes (1^3/2)^ or {2si/2)-n: and the particle (1/7/2)^- 

For the heavier isotopes, the situation is more complicated because the unpaired levels can be more than one, and 
because the levels involved in the gap calculation could be both of the same type (j> or j<). This is the case for Z 
or N=28 or 50, but not for Z — 40. For this reason, in the case of protons, the values of Ag are negative for Ni, Sn 
and Pb isotopes and positive for ^°Zr. 

For the neutrons, panel (a), the situation is more difficult to discuss since the nuclei we have investigated are not 
isotones. For each set of isotopes there are different s.p. levels related to the neutron gap, therefore the effect of the 
tensor term can be different for each nucleus considered. This is the reason of the oscillating behaviour we observe in 
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Figure 8; Differences between the energy gaps calculated with and without tensor force, for neutrons, panel (a), and protons, 
panel (b). The results obtained with the DIS type interactions are indicated by the open squares, and those obtained with the 
DIM type interactions by the solid circles. The arrows indicate those isotopes where the tensor effect is expected to be zero. 
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Figure 9: Single particle levels around the Fermi surface which change their order when the tensor force is used. 



the Ni isotopes. 

Since oxygen and calcium isotopes are spin saturated in protons, the results presented in the panel (a) of Fig. [5] 
for these nuclei, are produced by the interaction of an unpaired neutron s.p. level of j> type with the neutron s.p. 
levels just below and above the Fermi surface. The comparison of the results of these nuclei shown in the two panels 
indicates that the effect of the tensor interaction between like nucleons is smaller than that between neutrons and 
protons, reflecting the fact that we have used a tensor-isospin term in the interaction. 

We have seen that the tensor term affects more the s.p. than the bulk properties of the nucleus. The effects of the 
tensor force on the s.p. levels can modify their relative order. If this occurs for the levels near the Fermi surface, the 
spin of even-odd nuclei neighboring the nuclei investigated, which is determined by the last unpaired nucleon, should 
be modified. We found some cases where this happens and few of them are presented in Fig. ^ where we show the 
evolution of the states near the Fermi surface for the different interactions we have used. 

In the *^Ca case we consider the two proton states below the Fermi surface. The (lc?3/2)7r is a j< state, therefore 
its energy is lowered by the tensor, while that of the {2si/2)tt state remains essentially unchanged. This effect inverts 
the order of the two states, as we observe in the results of the DIST column. The calculations done with the DIM 
interaction give the two proton states in an inverted order with respect to that of the DIS one. In this case the tensor 
effect enlarges the energy difference between the two states. Since '*^Ca and ^®Ni are mirror nuclei, we expected an 
analogous effect for the two neutron hole states in ^^Ni. This effect is present in our calculations as it is shown in the 
figure. 

Since all the spin-orbit partner levels are occupied for both protons and neutrons, we do not expect tensor effects 
in ^°Ca. This is what we observe in Fig. [SI We show the results for this nucleus since the order of the proton hole 
states obtained with the DIS and DIST interactions is inverted with respect to that obtained with the DIM and 
DIMT forces. 

In the ^*Ni case, the tensor force produces a large lowering of the {1/5/2)7!- level, and this, for the DIM interaction, 
generates an inversion with the (2p3/2)7r level. The figure shows that the tensor force reduces the energy gap between 
(1/5/2)77 and (1/7/2)77 levels. 

For the ^'^''Sn nucleus we show both proton and neutron cases, since in both situations we obtain the inversion 
of the I57/2 and 2^5/2 levels, when the DIST interaction is used. The tensor force lowers the energies of the 137/2 
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levels, which are j< states, and enhances those of the 2^5/2 levels, which are j> states. We obtain a large effect on 
the energies of the I57/2 levels while the modification of the 2^5/2 energy is minimal. 

We observe an inversion of the order of the (2^3/2)1/ and (851/2)1/ levels in ^^''Sn, when the DIST interaction is 
used. The tensor lowers the energy of the (2^3/2)1^ state, of j< type, and has no effect on the (851/2)1/ state. The 
effect is present also when the DIMT interaction is used, but it is not large enough to invert the order of the states. 

In ^^^Sn we obtain an inversion of the order of the (l(77/2)7r and (2(i3/2)7r levels when both DIST and DIMT 
interactions are used. Also in this case we observe the tensor effect predicted by the Otsuka's scheme. The energy of 
the (137/2)77 state, of j< type, is lowered, while that of the (2^3/2)77 state, of j> type, is enhanced. 

To summarize the results presented in this section, we may say that the tensor effects have remarkable consequences 
on those observables which we relate to s.p. properties of the nuclear system, such as s.p. energies, gaps and spin of 
the system. 



IV. RPA CALCULATIONS 



In the previous section we have presented the results regarding the ground state properties of some spherical nuclei. 
In this section we show the results we have obtained for the excited states of these nuclei by doing RPA calculations. 
In this study, we have also considered ^^C, a well studied nucleus from both experimental and theoretical points of 
view, despite the fact that its ground state contains deformations [39i] . This feature is less important in the description 
of the excitation spectrum than for the ground state observables. 

The input required by any RPA calculations is composed by the s.p. basis and the effective interaction. The 
results shown in this section have been obtained by using the s.p. basis and the effective interaction provided by HE 
calculations. In the RPA jargon this procedure is called self-consistent to distinguish it from the phenomenological 
one described in Sec. ID To be precise, we do not strictly use the same interaction in HP and RPA, since in the latter 
case we neglect the Coulomb and spin-orbit terms of the interaction. We should keep in mind this fact, even though 
there are indications that, in RPA calculations, these two terms of the interaction produce very small effects that, in 
addition, have the tendency of canceling with each other (isj . 

The strategy of our investigation is analogous to that used in the previous section, i.e. we compare results obtained 
by using interactions with and without tensor terms. While for the HP calculations this comparison provides a clear 
indication of the effects of the tensor force, in the RPA case the situation is more complicated. In effect, the differences 
in the RPA results can be due to genuine effects of the tensor interaction in the RPA calculations, and also to the 
different values of the s.p. energies and wave functions provided by the HP results. 

We have disentangled the effects produced by these two different sources by presenting results where the tensor 
force is switched on and off in both HP and RPA calculations. To distinguish the results of the different type of 
calculations we have indicated with 1^^^ the RPA excitation energies, where the first subindex, a, refers to the 
interaction used in the HP calculation, and the second subindex, b, to that used in the RPA one. The a and b labels 
can be t if the interaction includes the tensor term, i.e. if the DIST or the DIMT interaction has been used, and n 
otherwise. Por example, ui^^^ indicates the excitation energy obtained with a HP calculation without tensor and a 
RPA calculation with the tensor force. We present also results obtained by switching off the residual interaction. We 
label these results as w^^^, where the superscript IPM means Independent Particle Model. The energies produced in 
this type of calculations need only one subindex which indicates the presence, t, or the absence, n, of the tensor force 
in the HP calculation. 

In our study we have investigated various types of multipole excitations, and, in agreement with the results of Refs. 
[Tsl [22l . [26j , we found that the natural parity excitations are practically insensitive to the effects of the tensor force. 
Por this reason, we present here only results obtained for the unnatural parity excitations. 

We start or discussion by presenting the results related to the excitation of the first 0~ state in the various isotopes 
we are studying. In the phenomenological calculations of Sec. [H] we have pointed out the large sensitivity of the 
excitation energy of these states to the tensor force. Por this reason we have chosen the energy of this excitation in 
^^O to select the strengths of our tensor forces. 

Por the nuclei where the 0~ excitation energies have been experimentally clearly identified p?! . [28| . we show in 
TableUthe excitation energies obtained in our RPA calculations. The tensor interaction lowers the excitation energies 
in all the cases. The perfect agreement of the ^^O case is obtained by construction, but we have improvements also in 
the *°Ca, ''^Ca and ^°®Pb cases. The only worsening produced by the inclusion of the tensor term is that of the ^^C 
nucleus, which we know to be a difficult nucleus to describe within our theoretical framework tailored to spherical 
systems. 

A general view of the tensor effects in all the nuclei we have considered is given in Pig. (TU] In panel (a) we show 
the differences wj^^-* — uj^^^ for the first 0~ states in all the nuclei under investigation. The open squares show the 
results for the interactions of DIS type, while the solid circles those of the DIM one. To separate the tensor effects in 
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Table I: Energies, in MeV, of the first excited state for those isotopes we have studied in this work where the energy values 
have been measured [53, [2^ . The theoretical energies are obtained by doing RPA calculations with different interactions. 



RPA calculations from those produced by the change in the s.p. configuration space, we show in panel (b) the energy 
differences ujI^^-oj^^^. 

In the RPA calculations, the results obtained with the tensor force are always lower than those obtained without 
it, and this produces negative values of the differences w^^^ — wj^jf'^. The only exception to this general trend is 
that of the ^^O nucleus in the case of the DIM interaction. In reality we observe in panel (b) that the differences 
^IPM _ (^IPM Jqj. ^j^jg nucleus are positive and larger than those found in the RPA calculation, indicating that also in 
this case the tensor force in the RPA calculation lowers the energy value of the first 0^ state. 

The results presented in the two panels do not show any correlation. This indicates that the effect shown in the panel 
(a) of the figure is a genuine effect of the tensor force on the RPA calculations, which, in the case of the 0~ excitation, is 
always attractive, even in the ^"^O case, confirming the results we have obtained with the phenomenological calculations 
and shown in Fig. [2l 

We have investigated the effects of the tensor force on multipole excitations with strong isoscalar (IS) or isovector 
(IV) character. To identify well these states we have considered only isotopes with N — Z. In these nuclei, we 
have selected those multipole excitations composed mainly by two identical particle-hole pairs in terms of angular 
momentum coupling, but one for protons and the other one for neutrons. This is the ideal situation to produce IS 
and IV partners levels. In the first case proton and neutron excitations are in phase, while in the second one are off 
phase. In our RPA results the IS and IV partner states are easily identifiable by observing the relative phases of the 
RPA amplitudes of the main particle-hole pairs. 

We have presented in Table |TT] the energies of the excited states we have investigated, and we compare them with 
the available experimental values [13] • All the results indicate that the energies of the IV excitations are larger than 
those of the IS ones, as it is experimentally well established. In Ref. [1^ we obtained for the calculations done with 
the DIS interaction an opposite behaviour. These results were wrong, since we found an error in the treatment of the 




Figure 10: Differences between the energies of the first 0~ excited state obtained with and without tensor forces, for the nuclei 
under investigation. DIS (DIM) results are plotted with open squares (solid circles). In the panels (a) and (b) we show the 
results obtained with RPA and IPM calculations, respectively. 
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Table II: Excitation energies, in MeV, in nuclei with Z = N, for different multipoles, where we have identified isoscalar (IS) 
and isovector (IV) character. The experimental values are taken from Ref . [i^l . 



exchange part of the density dependent terms of the DIS force. 

The effects of the tensor terms are better presented in Fig. [TT] In the panels (a) and (b) we show the energy 
differences uj^^^ — uj^^^ for IS (open squares) and IV (solid circles) states. The results found for the various 
calculations have rather similar behaviours. The tensor effects are smaller on the IV states. In the panels (c) and (d) 
we show the differences between the energies of the IV and IS states for each multipole we have considered. In these 
panels the open squares show the results obtained without tensor interaction, while the solid circles include the tensor 
in both HF and RPA calculations. The tensor force always increases these differences. In general, this enhancement 
is larger for DIS than for DIM and it is worth pointing out that the effect is relatively large for the three 1+ states 
studied (in ^^C and ^^Ni and ^"^Sn) in case of the DIS interaction. 

We have studied the effects of the tensor force on the electron scattering responses. A detailed presentation of 
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Figure 11: In the panels (a) and (b) we show the differences between the energies calculated with and without tensor terms for 
the IS and IV states listed in Table [III In the panels (c) and (d) we show the differences between the energies of the IV and IS 
states obtained in fully self-consistent RPA calculations with and without tensor forces. In the panels (a) and (c) we show the 
results obtained with the DIS type forces, and in the other two panels those obtained with the DIM type forces. 
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Figure 12: Inelastic electron scattering transverse response as a function of the effective momentum transfer [45|, calculated 
by using RP A wavefunctions obtained in fully self-consistent calculations done with different interactions. The data are from 
Refs. SES. 



these results would require a discussion for each specific excited state. We plan to make this discussion in the future. 
At the moment we can summarize the main and general feature we have observed by saying that the effects of the 
tensor force are larger on the IS excitations than on the IV ones. As example, we show in Fig. [T^] the inclusive 
inelastic electron scattering transverse responses [45| as a function of the effective momentum transfer, for the 1+ 
isospin excitation doublet in ^'^C, and the 2~ doublet in *°Ca, and we compare them with the available experimental 
data [i^, 113. It is evident that the effects of the tensor are greater on the IS than on the IV excited states. 

The explanation of these facts is related to the structure of the electromagnetic excitation of the unnatural parity 
states, which is dominated by the magnetization current [1^. The magnetization current depends on the anomalous 
magnetic moment of the nucleon, which has different sign for protons and neutrons. Since in the IS excitations 
the main proton and neutron RPA amplitudes have the same sign, the proton and neutron magnetization currents 
subtract with each other. In the IV excitations the effect is reversed. Small differences in the proton and neutron 
structure of the RPA wavefunctions are emphasized by the difference and hidden in the sum of the proton and neutron 
contributions. 

A remarkable sensitivity to the tensor force of the 1+ excitation of ''^Ca has been pointed out in Ref. |15|]. The two 
tensor forces used together with the Skyrme interactions in that article produce opposite results. We have investigated 
how our tensor forces affects the energy of the first 1"*" excitation in the nuclei with different number of protons and 
neutrons we have considered. We have summarized in Fig. [T2] our results. In this figure we show the differences 
between the 1+ excitation energies obtained in different type of calculations. The states selected are those with 
the largest _B-value and all of them are dominated by a neutron p-h configuration. Solid circles indicate the energy 
differences wj^^ — wj,^^. With the black squares we have shown the energy differences uj^^'^ — indicating the 

effect of the tensor forces in the RPA calculation. Finally, the open squares show the energy differences wj^^^ — Wj'^jf^, 
which are linked to the global effect of the tensor in our self-consistent calculations. 

The results shown in the figure indicate that the presence of the tensor force changes the s.p. energies in a way that 
the energy of the excitation is reduced. This is shown by the fact that all the energy differences have negative sign. 
These results are in agreement with the findings of Cao et al. [l^ for the T44 interaction [l^ and they have opposite 
behaviour of those found by the same authors for the modified SLy5 interaction [3] • In their calculations this is due 
to the change of the overall sign of the tensor term for the N-N pairs with the same isospin (tttt or i/iy pairs). 

As we have observed in all the results presented so far, the tensor effects are smaller for the DIM interaction. The 
effect of the tensor force is almost zero for the two nuclei where all the s.p. spin-orbit partner levels are occupied, 
i.e. ^^0 and ^"Ca. With the exception of these two cases, the differences ujft^^ — i^nn^ ^-^^ noticeably larger than 

, ,IPM , ,IPM 

To investigate the effect of the tensor force in RPA, we have calculated the energy differences w^^^ — 1^^^^ where 
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Figure 13: Differences between the energies of the 1^ excitations having the largest B(Ml)-value, in various nuclei. Solid 
circles indicate the results obtained in the IPM by using s.p. wave functions and energies obtained in HF calculations done 
with and without tensor. The solid squares show the results obtained in RPA calculations with and without tensor force, using 
the s.p. wave functions and energies obtained in HF with tensor force. The crosses show the sum of the two previous results. 
Open squares show the energy differences obtained by doing complete self-consistent calculations where the tensor force is used 
or not in both HF and RPA calculations. 



the same s.p. basis is used for the two RPA calculations. These results are presented in Fig. [12] by solid squares. We 
added incoherently the OJ^^ — oo^^^ and the w^^^ — Wtn^^ results and we obtained the results plotted with crosses 
which reproduce very well the results of the complete calculation cj^^^ — w^jf^. These results validate the assumption 
used by Cao et al. in Ref. [Tsj. 

We conclude this section by presenting the results obtained for the 1+ excitation of ^''^Pb which has attracted 
great attention in the past [isl - fsH . and whose interest has been renewed by recent measurements |52|. Energies and 
B-values for the first two 1+ excited states obtained with the four interactions are shown in Table IIIII We observe 
that the tensor lowers the values of the energies of both states. The B{M\) value of the first state is lowered by the 
tensor force, while that of the second excited state, which is the state included in Fig. [T31 is increased. 



exp PIS DIST DIM DIMT 

E{1+) [MeV] 5.85 7.80 4.76 6.50 4.82 

g(Ml)i 2.0 5.08 2.41 2.33 1.80 

E{lt) [MeV] 7.1-8.7 10.15 8.06 9.42 8.38 

5(M1)2 16.0 (17.9) 29.63 32.84 31.46 32.26 



Table III: Energies and B{M1) values, expressed in terms of nuclear magnetons, of the first two 1"*" excitations in Pb obtained 
by using different interactions. The experimental values are taken from Ref. [53 ]. 



The detailed analysis of the results shows that the effects of the tensor force on the energy values improve the 
agreement with the experimental values. The situation regarding the B{M\) values is more complex. The tensor 
force lowers the B{M1) values of the first state, and this improves the agreement with the experiment. The situation 
is reversed for the second excited state. Our results are compatible with those found in literature [13, IH, H^j • The 
large difference between the theoretical and experimental B{M\) values of the second state can be attributed to the 
limitations of the RPA which is unable to describe the large fragmentation of the second 1+ state. 

We show in Fig. [T4|the inelastic electron scattering responses [45|, calculated for these two states and we compare 
them with the data of Ref. [i^ . We observe that there is a larger sensitivity to the tensor force in the case of the 
first excited state. 
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Figure 14: Inelastic electron scattering transverse responses as a function of the effective momentum transfer [4^, calculated 
by using RPA wavefunctions obtained in fully self-consistent calculations done with different interactions. In panel (a) we show 
the results obtained with the RPA wavefunctions of the lowest excited states. The RPA wave functions used to obtain the 
results shown in panel (b) are those of the second excited state. The experimental data are from Ref. 49]. 



V. SUMMARY AND CONCLUSIONS 



We have studied the effects of the tensor terms of interactions used in HF and RPA theories. We considered only 
finite-range terms in the tensor-isospin channel, because they are the most important terms in microscopic interactions. 
We constructed these tensor terms by multiplying the analogous term of the microscopic Argonne V18 interaction by 
a function which modifies its behaviour at short internucleonic distances. This function contains a single parameter 
whose value determines the strength of the tensor force. 

To determine this strength, we searched for a global observable particularly sensitive to the tensor force, and, by 
using phenomenological RPA calculations, we found it in the excitation of the 0~ states. We have also observed 
that the tensor force strongly influences the values of the s.p. energies. We constructed new effective interactions by 
adding a tensor-isospin term to finite range interactions, and we selected the strength of the tensor and spin-orbit 
forces to reproduce, in the nucleus, the experimental energy of the first 0~ excited state, and the splitting between 
the s.p. energies of the neutron 1^3/2 and 1^3/2 levels. We based our work on the DIS and DIM parameterizations 
of the Gogny interaction, and we called DIST and DIMT the two new effective interactions constructed with a 
recursive procedure where HF and RPA calculations have been repeated until both experimental quantities have been 
reproduced. 

The study of the effects of these tensor forces has been done by comparing the results obtained in HF and RPA 
calculations, using interactions with and without tensor terms. We have repeated each calculation with both DIS 
and DIM type of effective interactions to extract genuine tensor effects from those related to the peculiarities of a 
specific choice of the force parameters. We always found great similarities between the results obtained with the two 
different type of parameterizations. Since the strength of the tensor force in DIMT is weaker than that of DIST, we 
found smaller tensor effects in the results obtained with DIMT than in those for DIST. We have done calculations 
for a set of spherical nuclei chosen such as all the s.p. levels below the Fermi surface are fully occupied. With this 
choice we avoided the effects of the deformation and we minimized those of the pairing. 

Our HF calculations indicate that tensor forces do not produce sensitive effects on bulk observables such as binding 
energies, radii and density distributions. The effects of the tensor force on quantities related to the s.p. properties 
of our theoretical approach are certainly more remarkable. We have calculated energy splitting between spin-orbit 
partner levels, energy gaps between the s.p. states just above and below the Fermi surface, and we found noticeable 
effects produced by the tensor force. Also the ordering of the s.p. levels around the Fermi surface, which determines 
the spin of the even odd nuclei neighboring those we have studied, is strongly infiuenced by the presence of the tensor 
force in some cases. We could explain all our results within the picture proposed by Otsuka and collaborators ^9i.[l0|. 
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eventually by extending it to consider the interaction between nucleons of the same type. 

We have already mentioned in the introduction that the clear identification of the tensor effects on the s.p. properties 
from the comparison with experiment is rather problematic. In the MF picture the effects of the tensor force are of the 
same size of those produced by the spin-orbit force. Beyond the MF model, it is well known that s.p. properties are 
strongly affected by correlation effects, the most important ones, for the quantities of our interest, are those related 
to the coupling of s.p. states with low- lying collective vibrations jssj ]. 

The study of tensor effects in the excitation spectra is more complex, since the effects of the tensor force in the RPA 
calculations add to those already present in the HF calculation that produces the s.p. bases. In our investigation we 
have disentangled the effects coming from these two different sources. 

We have verified the well know fact [13, that natural parity excitations are essentially unaffected by the tensor 
force. For this reason we have presented results regarding unnatural parity excitations only. We started our investi- 
gation by studying the 0~ excitations in various nuclei, and we found that the excitation energies obtained with the 
tensor forces are always smaller than those obtained without it. 

We studied the different role played by the tensor force in IS and IV type of excitations. To identify clearly these 
different type of excitation modes, we have considered nuclei with equal number of protons and neutrons, and states 
dominated by particle-hole transitions with the same angular momentum coupling for both protons and neutrons. 
The energies of the IV modes are always greater than those of the IS modes with or without tensor. These results 
agree with the experimental observations. We found that the IS excitations are more sensitive to the tensor force 
than the IV ones. The tensor force increases the energy difference between IS and IV excitations. 

We made a systematic study of the effects of the tensor force on the excitation of the 1+ states in nuclei with 
different number of protons and neutrons. We have considered the states showing the largest _B(Ml)-values, and we 
found that the tensor force consistently lowers the values of their excitation energies. This is essentially obtained as 
the incoherent sum of the effect generated by the HF calculation and that obtained by the RPA. 

The study of the excitation of the first two 1+ excited states in ^°*Pb indicates that the lower energy state is more 
sensitive to the tensor force than the other one. The presence of the tensor force modifies energies and i?(Afl)- values, 
and slightly improves the agreement with the experimental data [s^ , even though the limitations of the RPA theory 
do not allow a description of the fragmentation of the strength of the higher energy state. 

We have restricted our study to charge conserving nuclear excitations. There are indications that the effects of the 
tensor-isospin force are more relevant on the charge-exchange excitations (53 - [56| . The presence of the tensor forces 
seems to be relevant also in the description of the reactions between heavy nuclei [STj . 

We think that the accuracy required today by self-consistent effective theories requires the use of interactions 
containing tensor terms. 
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